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It is shown that almost aU known solutions of the kind mentioned in the title are easily derived 
in a unified manner when a simple ansatz is imposed on the metric. The Whittaker solution is an 
exception, replaced by a new solution with the same equation of state. 
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I. INTRODUCTION 

T-H . 

f^ ' Einstein's equations for static and spherically symmetric perfect fluids have been studied by many authors [1-3]. 

^^ I The unknown functions are more than the equations, therefore, either an equation of state should be prescribed or 
Cn . some simplifying ansatz assumed. Schwarzschild found the flrst solution [4] in 1916. It was a degenerate case (n = 0) 
of the linear equation of state 



o 



X 



np + po (1) 



" , where p and p are the density and the pressure of the fluid and n, po are constant parameters. A year later, two other 

^^ ■ solutions satisfying Eq. (1) were found, describing static cosmological models [5,6]. Tolman [7] used eight simplifying 

►^ , ansatze for the metric components to find new solutions. Some limiting cases of these satisfy Eq. (1) too, leading 

fvj ' to a simple, but singular solution for arbitrary n and po = [3]. Klein rediscovered this solution first for the case 

(T^ [ of incoherent radiation n = 3 [8] and then for general n [9] by studying systematically the equation of state (1) 

^^ ■ with vanishing po (the so-called 7-law). This Klein- Tolman (KT) solution was rediscovered later three times more 

C ■ [10-12]. Meanwhile, two new solutions with pq j^ appeared in 1968, developing further the idea of Schwarzschild's 

^^ incompressible sphere. Whittaker argued that the gravitational mass density p + 3p should be constant [13], while 

^«. ' Buchdahl and Land demanded that the speed of sound should be the largest possible, namely, equal to the speed of 

^ ■ light [14]. 

^ I The main stream of papers, however, utilized simplifying assumptions, which often lead to complicated equations 
of state, if any. Already Tolman used the ansatz 



I 

W). e-^ = a + br^ (2) 



where a and b are constants, < a < 1, and e^ is the g^r component of the metric. A thorough study of this ansatz 
was performed by Kuchowicz [15,16] who found seven different cases in total, many of them having rather involved 
expressions for the pressure and the density. He didn't study whether his solutions satisfy Eq. (1). 

What happens when we impose on the field equations both Eqs. (1) and (2)? An overdetermined system of 
differential equations is obtained which may be devoid of any non-trivial solutions. The real answer, derived in this 
paper, is somewhat surprising. All known explicit solutions with a linear equation of state satisfy also Eq. (2), except 
for the Whittaker solution. We find an additional solution with n — —3 instead, which appears to be new. 

It is much easier to impose both requirements directly on the field equations, than to study case by case the results 
of Ref. [16]. In this way we derive in a unified and simple manner almost all known solutions with linear equation of 
state, including a new hybrid between the Einstein static universe (ESU) and the KT solution. This programme is 
executed in the following section. The last section contains some discussion. 

II. FIELD EQUATIONS AND THEIR SOLUTIONS 

The metric of a static spherically symmetric spacetime reads 

ds^ = e^dt^ - e^dr^ - r^dn"^ (3) 

where dO^ is the metric on the two-sphere and v, A depend on r. The system of equations is 



p=^e- + l(l-e-) (4) 

p' + ^{p + p)„' = (6) 

where the prime means a r-derivative and one of the field equations has been replaced by the contracted Bianchi 
identity (6). Let us impose Eq. (2). Then Eq. (4) gives a simple expression for the density 

1 — a 

P=-'^b+-^ (7) 

which is regular at the centre if a = 1. When Eq. (1) holds, Eq. (6) may be integrated. Two different cases arise, 
n ~ —I and n 7^ — 1. 

Let us discuss first the case n = — 1. Then Eq. (6) yields 

p = --poi^ + h (8) 

with h being an integration constant. Plugging Eq. (8) into Eq. (1) and the resulting expression for p into Eq. (7) 
gives 

1 1 — a , , 

-poiy = ft - po - 36 H 5— (9) 

There are two possibilities. If po = then from Eq. (8) p = h, from Eq. (1) p = —p and from Eq. (9) a = 1, 6 = h/S. 
Therefore 

e-'^l + y' (10) 

while Eq. (5) gives v after a simple integration: v ~ —X. This is the de Sitter solution. 
If po 7^ then the insertion of Eqs. (2,8,9) into Eq. (5) yields the relation 

po (po + 26) r^ = -4 (1 - a) (a + 6r2) (11) 

Obviously, a = 1 and 6 = — po/2. From Eq. (7) p = 3po/2 and from Eq. (1) p = — po/2. Eq. (8) gives i' = 1 + 2/i/po = 
const, while Eq. (2) becomes 

e'^ = l- Plr^ (12) 

2 ^ ' 

The pressure and the density satisfy, in fact, the equation of state p + 3p = 0. Hence, we obtain ESU, although 
initially p = —p + po . 

Let us continue with the generic case n 7^ — 1. Now the integration of Eq. (6) brings a different answer 

p = Pie-'^'^-po (13) 

where pi is a constant of integration, while pg = Po/ (^ + 1)- Proceeding like we have done with Eq. (8) the analog 
of Eq. (9) is obtained 

n+i 1 — a 
npie 2 -npo + Po = -36H — (14) 

Again, there are two possibilities. If n = we have a = 1, 6 = — po/3 and from Eq. (1) p = po = const. The metric 
component v is obtained not from Eq. (14) but from Eq. (5), the last equation to be satisfied. For a general n it 
becomes a linear equation for y = e"^-'^ : 



In our case it reduces to 



and is easily integrated 



—r (a + 6r^) y' = [l - a - (po + b) r^] y + pir"^ 



2 (^1 - y?" J y = -^rpoy+pir 



(15) 



(16) 



2po 



(17) 



e-A = 1 _ ^,2 



(18) 



which, after some change of notation, coincides with solution III of Tolnian. This is the Schwarzschild interior solution. 
If n 7^ both Eqs. (14) and (15) determine ly and we must study their compatibility. Eq. (14) may be written as 



- ^A + — 
y r 



(19) 



A = -- 



-FTT)^t^° + '^" + '^'^ 



(20) 



B = 



1-a 
npi 



(21) 



When these relations are put into Eq. (15), a polynomial of second degree in r^ arises, which must vanish identically. 
Three compatibility conditions result 



{I -a) 



4n+ (n+ 1)' 
{n + lf 



a-1 



= 



(22) 



A [po + {n + 3)b]=0 



(23) 



B [in + 1) Po + 4nfe + (n + 1) (n + 3) b] - (n + if {1 - a) A ^ 



(24) 



Eq. (22) has two solutions for a. Let a — I. Then B = 0, j^ is a constant and p = —36. Eq. (24) is also satisfied, 
while Eq. (23) provides two subcases. If ^ = then from Eq. (19) y = oo which is unacceptable. li po ~ — {n + 3) b 
then Eq. (20) shows that A ^ when b ^ 0. Eq. (1) gives p = b, hence the equation of state is p + 3p = for any n. 
We again obtain the ESU with constant i' and 



e"^ = l + br^ 

This is the same as Eq. (12), since both can be written as e^^ = 1 +pr'^. 
The last case which remains is 

{n+lf 



4n + (n + 1) 



(25) 



(26) 



where a ^ 1 since n 7^ 0. Eq. (23), like before, offers two choices. 

a) A = (po = —3 (n + 1) b). Then Eq. (24) simplifies to (n — 1) 5 = 0. There are two subcases: 

al) n = 1, 6 7^ 0. Then po = -66, Eq. (26) gives a = 1/2, Eq. (21) gives B = l/2pi and Eqs. (2,19,7,1) give 

respectively 



-A 1 /i Po 2 
e = - 1 —r 



(27) 



e'' = 2pir'^ 



P-2[P' + ^ 



P = P + Pa 



This is precisely the solution of Buchdahl and Land [14] after e^ is rescaled. 
a2) h = {). Then po = and we easily find 

{n+lf 



(28) 
(29) 
(30) 

(31) 
(32) 

(33) 



This is the KT solution in curvature coordinates [3,8-12]. Eq. (31) shows that the solution exists when 4n+ (n + 1) > 
0. This condition is fulfilled for n outside the interval (-5.83 = -3 - 2\/2, -3 + 2^2 = -0.17). When n < -5.83 or 
—0.17 < n < the solution has p and p of different signs. The physically realistic range is 1 < n < oo, which includes 
the important cases of stiff fluid (n = 1) and incoherent radiation (n = 3). The KT solution is singular at r = but 
it can serve as an interior solution for r > tq > in multi-layered perfect fluid models. It also represents a focal point 
for a regular solution, which has been studied numerically [8,17,18, 19]. 

h) po = — {n + 3) b. Since a ^ 1, Eq. (24) reduces to {n + '3)b = 0. The two subcases are 

bl) 6 = 0. Then ^ = and we go back to case a). 

b2) n = -3. Then po = 0, a = -1/2, p + 3p = and 
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4n+(n+l)^ 
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l^br^ 



(34) 



Pie 



1 

2^ 



(35) 



P 



-36 



2r2 



(36) 



The constant pi may be set to one by a time rescaling. This appears to be a new solution - a hybrid between ESU 
with its constant pressure and density and the KT solution. Its pressure, density and metric component gtt are linear 
combinations of those two solutions, although Eqs. (4-6) form a highly non-linear system. When 6 = we obtain the 
KT solution for n — —3 which, in fact, does not exist. When b y^ 0, e^ and e'^ are negative for small r. However, if we 



choose 



1/26 and 6 > 0, then the solution exists for r > tq, has positive pressure, negative density and represents 



a regular companion of ESU. 

Whittaker has solved the general case n = —3 when po may be different from zero and the ansatz (2) is not imposed 
[13]. However, the above solution has been missed. The reason is the following. Eqs. (1,4,5) lead to the relation 



(6^+^^)' = A:re2(^+'') 



where k is a constant. The integration of this equation yields 



c — kr^ 



(37) 



(38) 



c being a constant. Whittaker assumes that k = ca, where a is yet another constant. This choice excludes the above 
solution which has c = 0, fc = — 2pi ^ and in principle may be regained as a limiting case. However, when one sets 
Po = in the results of Ref. [13] only ESU is obtained. It should be noticed that the singular term in Eq. (35) is 
different from the universal Schwarcschild term m/r which was correctly discarded from the interior solution. 



III. DISCUSSION 

In this paper wc have studied the gravitation of static spherically symmetric perfect fluid solutions with a linear 
equation of state. We imposed the additional requirement (2), making the system of equations overdetermined. 
Amazingly, it gives six non-trivial solutions: the Schwarzschild interior solution for an incompressible sphere, the de 
Sitter eosmological solution, the Einstein static universe, the stiff fluid solution of Buchdahl and Land, the Klein- 
Tolman singular solution and a new solution with p + 3p = 0. Among the known solutions with linear equation of state 
only the Whittaker solution does not satisfy Eq. (2). The reason is that it satisfies Eq. (38) and if Eq (2) is required 
too, ly has the only option to remain a constant, which is too restrictive. These solutions have been found originally 
starting from a variety of physical considerations and following quite different mathematical approaches. Here they 
are derived in a unified way as a result of straightforward calculations and a new solution is added. Roughly speaking, 
it is a linear combination between ESU and the KT solution, which is unusual for a non-linear system of differential 
equations. It closes the list of solutions satisfying both Eqs. (1) and (2). The equation of state p + 3p = const is 
unrealistic because the speed of sound is negative. However, models which satisfy it have interesting properties, are 
often integrable and occupy important places in different classification schemes. This concerns one whole class of 
algebraically special perfect fluids of Petrov type II or D [1] (Eq.(29.20)), a rigidly rotating solution of type D [20], 
the well-known Wahlquist solution [21-22] which has the Whittaker solution as a static limit and the Kerr metric as 
a vacuum subcase, and the recently found Wahlquist-Newman solution [23] which reduces to a number of well-known 
solutions in appropriate limits. An interesting phenomenon is the effective change of the equation of state by which 
ESU is obtained first from p = — p + po and then from p = np + pQ. When pa = a one-parameter sequence of 
KT solutions exists which serve as focal points for regular numerical solutions. When po 7^ only two discrete cases 
satisfy Eq. (2) for n — and n — I, the third known one (n = —3) drops out of this scheme. Presumably when one 
replaces the ansatz (2) by a different one, other integrable cases will appear, although the system is overdetermined. 
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